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ABSTRACT 

We examine deep XMM-Newton Reflection Grating Spectrometer (RGS) spectra from the 
cores of three X-ray bright cool core galaxy clusters, Abell 262, Abell 3581 and HCG 62. 
Each of the RGS spectra show Fe XVII emission lines indicating the presence of gas around 
0.5 keV. There is no evidence for O Vll emission which would imply gas at still cooler tem- 
peratures. The range in detected gas temperature in these objects is a factor of 3.7, 5.6 and 
2 for Abell 262, Abell 3581 and HCG 62, respectively. The coolest detected gas only has a 
volume filling fraction of 6 and 3 per cent for Abell 262 and Abell 3581, but is likely to be 
volume filling in HCG 62. Chandra spatially resolved spectroscopy confirms the low volume 
filling fractions of the cool gas in Abell 262 and Abell 3581, indicating this cool gas exists as 
cold blobs. Any volume heating mechanism aiming to prevent cooling would overheat the sur- 
roundings of the cool gas by a factor of 4. If the gas is radiatively cooling below 0.5 keV, it is 
cooling at a rate at least an order of magnitude below that at higher temperatures in Abell 262 
and Abell 3581 and two-orders of magnitude lower in HCG 62. The gas may be cooling 
non-radiatively through mixing in these cool blobs, where the energy released by cooling is 
emitted in the infrared. We find very good agreement between smooth particle inference mod- 
elling of the cluster and conventional spectral fitting. Comparing the temperature distribution 
from this analysis with that expected in a cooling flow, there appears to be a even larger break 
below 0.5 keV as compared with previous empirical descriptions of the deviations of cooling 
flow models. 
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1 INTRODUCTION 

The hot (few lO^K) intracluster medium (ICM) in galaxy clusters 
is primarily seen in emission in the X-ray waveband. The ICM con- 
tains most of the baryonic cluster mass. In a large fraction of clus- 
ters of galaxies, observations show that the X-ray surface bright- 
ness steeply rises towards their centres (e.g. [Stewart et al.||1984| ). 
In addition, spectrally-derived temperature profiles of clusters typ- 
ically show drops in temperatures by a factor of 2 or 3 from the 
outskirts (e.g. [Allen et al.|2001| ). As the X-ray surface brightness 
is proportional to the integral of the density squared along a line 
of sight, peaked surface brightness profiles imply that the gas in the 
cores is cooling much more rapidly than in the outskirts. In the cen- 
tre, the mean radiative cooling time, often calculated as the ratio of 
the luminosity of a region to its enthalpy, often drops below 1 Gyr. 

Assuming steady-state and the absence of heating, there will 
be material with short cooling times in the cluster centre cooling 
out of the X-ray emitting band. As the volume of this material be- 
comes much smaller, there must be a flow of material in to replace 
it to preserve the steady state. The steep surface brightness profiles 
imply rates of 10s to 1000s of solar masses per year cooling out of 
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the X-ray band (see |Fabian| 1 994| for a review). This cooling mate- 
rial would be expected to eventually give rise to star formation. 

This picture changed when high spectral resolution X-ray 
studies of nearby clusters of galaxies using the Reflection Grating 
Spectrometer (RGS) instruments on XMM-Newton revealed a lack 
of cool X-ray emitting gas in these objects ( Tamura et al. "200 1 a|b [ 
Peterson et a l. 2001 , Kaastra et al. 2001 , Sakelliou et al. 2002 ;1Pe^ 
terson et aL||2003| ). The main spectral indicator missing in these 
spectra are the emission lines of Fe XVII, which are strong from 
material between 0.15 and 0.8 keV (1.7 to 9.4 MK). 

The general accepted picture is that there is a lack of material 
below a factor 2 or 3 of the outer temperature, which matches the 
results from Chandra spatially -resolved temperature profiles. This 
is not completely correct, however. Deep observations of nearby 
clusters show a much larger range of temperature. Fe XVII emis- 
sion lines have been seen in Centaurus ( [Sanders et al.|2008'] ), show- 
ing a temperature range of more than 10, Abell 2204 ([Sanders 
et al.|2009b| ), showing a range of 15, 2A 0335-^096 ( [Sanders et al. 
2009a) ), showing a range of 8. Fe XVII emission was also detected 
in Abefl 262, M87 and NGC 533 in a large study by jPeterson et al.[ 
( [2003 1 ). There is much less material than would be expected to be 
seen in the case of steady- state radiative cooling without any heat- 
ing, however. 
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There is evidence that the central active galactic nuclei (AGN) 
in these objects can energetically prevent much of the cooling (see 
reviews by [Peterson & Fabian] |2006| and jMcNamara & Nulsen| 
[200 7"). This heating may be by the inflation of cavities (radio lobes) 
by the AGN or the subsequent dissipation of sound waves gener- 
ated by the inflation. Such cavities are found in almost all nearby 
clusters which require heating fP unn & Fabian|2006| . 

If AGNs are responsible for preventing cooling in cluster 
cores, there are still remaining issues. In Centaurus the cooli ng time 



of the lowest detected component is only 10^ yr (Sanders et al. 



[2008] ). If no cooling is taking place, then feedback must be able to 
operate on these timescales. Feedback must also be able to operate 
over the much longer cluster lifetime. Star formation in Centaurus 
must have occured slowly over the last 8 Gyr, or all must have been 
at earlier times ( [Sanders & Fabian|2006| l. 

We assume a Hubble constant of 70kms~^ Mpc~^ and use 
the relative Solar metallicities of [Anders & Grevesse| ( [1989l ). 



2 THE SAMPLE 

We are investigating cool gas in a sample of three nearby ob- 
jects with a range of mass scales. Table [T] shows the objects, red- 
shifts. Galactic absorptions, XMM-Newton observations, observa- 
tion dates and exposures. 

The Chandra data of Abell 262 were examined by [Blan-[ 
[ton et al.| ( |2004| and [Clarke et a l. ( 2009). The cluster lies in the 
Perseus supercluster and has a bolometric X-ray luminosity of 
4.4 X 10^^ erg s~^ from Einstein ( David et al. 19931 ). Its cD galaxy 
is NGC 708. There is extensive interaction between the radio source 
and X-ray emitting gas. There is a tunnel like region of lower X- 
ray surface brightness to the south west which is coincident with 
an extension of the radio source. Towards the east there are three 
clumps of radio emission. The innermost one is coincident with an 
X-ray cavity. There is an additional cavity to the east. The bright 
X-ray structures are in the same regions as the optical line emis- 
sion. JCT^eetari (2009) put a timescale on the outburst repetition 
timescale of the central nucleus of 28 Myr. They claim that the to- 
tal AGN output is sufficient to offset radiative cooling over several 
episodes of outbursts. The coolest X-ray gas seen by Chandra is 
around 0.8 keV ( Blanton et al.,2004| ). 

The short Chandra observation of the galaxy cluster 
Abell 3581, which contains the radio source PKS 1404-267, was 
examined by [Johnstone et al.[ ( [2005| ). They found that the X-ray 
temperature declines to around 0.4 of the maximum within the 
central 5 kpc and where the entropy drops to below lOkeVcm^. 
There are clear X-ray cavities located at the position of the radio 
source at 1.4 GHz. [White e"taL| ( [T997] ) give an X-ray luminosity of 
4.2 X 10^^ erg s~^ for the cluster, correcting to our cosmology. 

jMorita et al.[ ( 2006) examined Chandra and XMM-Newton 
data of the bright Hickson compact galaxy group HCG 62. They 
find two cavities in the X-ray emission. Two temperature compo- 
nents were required within the inner 2 arcmin (0.7 and 1.4 keV). 
The bolometric X-ray luminosity of the cluster is 5.6 x 10^^ erg s~^ 
( [Ponman & Bertram] 1993] ). 



reprocessed to apply the latest gain and charge transfer inefficiency 
correction and then filtered to remove bad grades. The improved 
background screening provided by VFAINT mode was applied to 
the observations of Abell 262 and Abell 3581. The background 
light curves of the level 2 event files were then filtered using the 
LC_CLEAN scripj^ provided by M. Markevitch to remove periods 
affected by flares. The final cleaned exposure times are shown in 
Table J2] Blank- sky background data sets available from the CXC 
were processed in the same way as the cluster observations and 
normalized to the source count rate in the 9 — 12 keV energy band. 

Shown in Fig. [T] are Chandra images of each of the three ob- 
jects. 



3.1 Deprojected quantities 

For each cluster, projected spectra were extracted in a series of con- 
centric annuli centred on the emission peak and then deprojected 
with DSDEPROJ ( [Sanders & Fabian[[2057l [Russell et ar[[2058] ) to 
remove the contributions from the outer cluster layers. We assume 
that the clusters are spherically symmetric in this analysis. The ra- 
dial bins were chosen to ensure a minimum of 3 000 deprojected 
counts in each deprojected spectrum. Point sources were identified 
using the CIAO algorithm WAVDETECT, visually confirmed and ex- 
cluded from the analysis. All spectra were analysed in the energy 
range 0.5 — 7keV and grouped with a minimum of 50 counts per 
spectral bin. Response and ancillary response files were generated 
for each cluster spectrum, weighted according to the number of 
counts between 0.5 and 7keV. The deprojected spectra were fit- 
ted in XSPEC version 12.5.0 ( [Arnau d 1996) with an absorbed ther- 
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19861 [Kaastra[[1992| 


Liedahl et al.[[1995|). The absorbing column 



3 CHANDRA DATA 

Chandra observations of the clusters (listed in Table |2]) were anal- 
ysed using CIAO version 4.1.2 with CALDB version 4.1.2 provided 
by the Chandra X-ray Center (CXC). The level 1 event files were 



density was fixed to the Galactic values given by I Kalberla et al.[ 
( 2005), shown in Table [T] The temperature, metallicity and model 
normalization were allowed to vary. We show deprojected temper- 
ature and electron density profiles in Fig.|2] 

The lines shown on the temperature and density panels in 
Fig. )2] are fr om a surface brightness deprojectio n (following [Fabian[ 
[et aljl981| ) assuming an NFW potential (Navarro et al.|1996|. We 
use a Metropolis-Hastings Markov Chain Monte Carlo (MCMC) 
method to iterate over the NFW parameters (concentration and 
r2oo) and the outer pressure in the deprojection. The MCMC cal- 
culates the of the NFW-predicted gas temperature profile to 
the observed deprojected spectral temperature profile to determine 
whether to move to the new position in the chain. The parameters 
are stepped individually in random order in the MCMC using a 
Gaussian proposal probability distribution. The width of the pro- 
posal distribution is dynamically adjusted to achieve a repeat frac- 
tion of 0.75 during the burn in period of 10 000 iterations. The lines 
show the median density and temperature profiles obtained with the 
MCMC. The input surface brightness profiles were created using 
radial bins a factor of three times smaller than the spectral bins for 
Abell 262 and HCG 62, and a factor of five for Abell 3581. 

From the set of potentials, pressures, densities and tempera- 
tures obtained from the surface brightness deprojection, we can ob- 
tain mean radiative cooling time profiles and cooling flow mass de- 
position rates. These are shown in the bottom two panels of Fig.j2] 
The median values from the MCMC are shown with the l-cr errors 
calculated from the 15.9 and 84.2 percentiles. The mass deposition 
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Table 1. Details of the targets and individual XMM -Newton observations. The exposure times given are for the RGSl instrument after cleaning. The absorption 
quoted is the Galactic absorption from^ Kalberla et al.]j2005) . References in superscript for redshifts are (l) |Struble & Rood|jl999) , (2) | Johnstone et al.|jl998t 
and (3) [Zabludoff & Mulchaey ( 2000 >. 



Cluster Redshift Absorption Observation Date Exposure 

(10^0 cm-2) (ks) 



Abell 262 


O.OI63I 


5.67 


0109980101 


2001-01-16 


26.2 








0504780101 


2007-07-12 


121.2 








0504780201 


2007-07-18 


41.4 








total 




188.8 


Abell 3581 


0.02 18^ 


4.36 


0205990101 


2004-01-29 


43.5 








0504780301 


2007-08-01 


116.3 








0504780401 


2007-08-03 


27.7 








total 




187.5 


HCG 62 


0.014533 


3.31 


0112270701 


2003-01-15 


12.4 








0504780501 


2007-06-26 


110.1 








0504780601 


2007-06-29 


33.6 








total 




156.6 




Figure 1. Chandra images of the cores of the three objects examined in this paper. Images have been smoothed with a 1 arcsec Gaussian. These images show 
the Chandra observations in Table |2] 



Table 2. Chandra datasets analysed. mass as a function of temperature or cooling time. These are plotted 
in Fig. [3] 



Target 


Obs-ID 


Date 


Clean exposure (ks) 


Abell 262 


7921 


2006-11-20 


110.7 


Abell 3581 


1650 


2001-06-07 


7.2 


HCG 62 


921 


2000-01-25 


48.4 



rates assume steady state and no heating and account for the gravi- 
tational contribution in the cooling. 

Shown in Table |3] are the best fitting NFW parameters and 
outer pressure in the deprojection. We also show the cooling radius, 
where the cooling time is approximately the time since z= I. The 
listed mass deposition rate is the cumulative mass deposition rate 
inside this radius. 

Note that the best fitting NFW concentrations are very high. 
This is likely to be because the Chandra observations concentrate 
on the region where the potential of the central cluster galaxy is 
most important. Only a very small cluster volume is being sampled. 
We also tried using a nonsingular isothermal model with the modi- 
fied Hubble law density approximation. The resulting density, tem- 
perature, cooling time and mass deposition rate profiles are almost 
identical with those computed using the NFW mass model. The me- 
dian velocity dispersion, a and King radius, tq, for the isothermal 
mass model are also shown in Table [3] 

Using these profiles we can calculate the distribution of gas 



3.2 Temperature maps 

We have calculated projected emission- weighted temperature maps 
of each of the clusters from the Chandra data (Fig.|4]). Regions were 
selected to have a signal-to-noise ratio of 22 (corresponding to 480 
counts in the absence of background) using the Contour Binning 
algorithm ( [Sanders |2006| ). The regions were chosen to have a max- 
imum ratio of the 2 between the longest and shortest dimensions. 
The spectra from the regions were fit between 0.5 and 7 keV with 
a single component APEC spectral model { [Smith et al.|2001| ) with 
free temperature, metallicity and normalization, minimizing the C- 
statistic in the fits. The Galactic absorption and redshift were fixed 
to the values given in Table [T] 

The temperature maps show a significant amount of structure. 
One interesting feature in Abell 262 is a 20 arcsec (6.7 kpc) radius 
ring of cool gas 3 arcmin (60 kpc) to the north east of the cluster 
centre. It is marked by arrows in the top panels of Fig. [4] The ring 
appears that it may be connected to the centre of the cluster with 
a cool ridge and the cool gas may extend northwards beyond the 
ring. The ring is not associated with any galaxies. 
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Radius (kpc) Radius (kpc) 

Figure 2. (Top left panel) Deprojected densities measured from Chandra spectral fitting (shown as points) and from surface brightness deprojection (shown as 
solid lines). (Top right panel) Deprojected temperature measured from Chandra spectral fitting (shown as points) and surface brightness deprojection (shown 
as lines). (Bottom left panel) Cooling time profiles measured from surface brightness deprojection (shown as small points) and RGS (shown as large points). 
(Bottom right panel) Cumulative mass deposition rates derived from Chandra surface brightness deprojection. 



Table 3. Surface brightness deprojection results. The values shown are the median from the NFW MCMC analysis. The error bars were calculated from the 
15.9 and 84.2 percentiles and are equivalent to l-cr errors. The values of a and ro were instead calculated using a nonsingular isothermal model (with the 
modified Hubble law form) in the MCMC analysis. 



Cluster 


Concentration 


^200 (Mpc) 


Pouter (logio erg cm 3) 


rcool (kpc) 


M(Mq yr-i) 


O" ( km s 1 ) 


ro (kpc) 


Abell 262 
Abell 3581 
HCG 62 


15.0^[9 
45.0+3-1 


342+0-012 
57+0-06 
0.271 ±0.008 


-10.74±0.01 
-11. 04 ±0.04 
-11.19±0.02 


39 
57 
36 


18±1 
49±6 
6.7±0.2 


256 ±3 
323 ±7 
222 ±3 


4.7±0.3 

13.1±1.6 

3.1±0.2 



4 RGS DATA ANALYSIS 

We extracted spectra from each observation listed in Table [T] us- 
ing RGSPROC, part of XMMSAS (version 8.0.0). For the main part 
of the analysis we extracted spectra from 95 per cent of the cross- 
dispersion point- spread-function (PSF) and within 90 per cent of 
the pulse-height distribution. The PSF extraction region corre- 
sponds to a roughly 100 arcsec strip across the cluster centre. 
As the cluster emission fills the RGS field of view, we generated 
template background spectra for each observation with RGSBKG- 
MODEL. Each of the observations for a target were processed using 
'attstyle' option in RGSPROC set to 'user' so that the spectra could 
be added together using RGSCOMBINE. As the wavelengths in the 
spectrum are dependent on the correct position of the source, we 
examined the Chandra observations by eye to find the coordinates 
of the peak diffuse X-ray emission. 

We applied time filters to remove the strongest flares. Al- 
though there are periods remaining in the observations containing 



relatively mild flares, the template background generator is able to 
generate appropriate backgrounds by adding together backgrounds 
with similar count off-axis count rates. After background subtrac- 
tion, each observation of each cluster matches the other observa- 
tions well at long wavelengths, where background is most impor- 
tant. We use a 7 to 26A wavelength range to minimize background 
effects. 

The observations for each cluster were combined with the 
RGSCOMBINE tool. The template backgrounds were combined with 
our own program which averaged them together, weighting by the 
exposure time of the respective foreground spectrum. 

Combined fluxed spectra from each of the targets, produced by 
RGSFLUXER and removing model backgrounds are shown in Fig. [5] 
In the lower panel we show similarly broadened model spectra at 
four different temperatures for comparison. Although the fluxed 
spectra cannot be used for spectral analysis, they show that the 
spectra contain Fe XVII emission lines indicating material around 
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Figure 3. Distribution of gas mas as a function of mean radiative cooling 
time and temperature. The distributions are calculated from the Chandra 
deprojection analysis (Fig.|2j. The points show values calculated from the 
RGS emission measures of the 0.544 keV component, assuming that they 
are in pressure equilibrium with the gas from the Chandra deprojection. 



0.7 keV (the ionisation fraction of Fe in Fe XVII peaks at 0.35 keV; 
[Mazzotta et al.|p"998| ). They do not, however, show evidence for 
O VII lines, indicative of gas below 0.25 keV (O Vll ion fraction 
peaks at 0.07 keV). 

The spectral features from the innermost region can be seen 
more clearly from spectra extracted using a smaller extraction re- 
gion. Shown in Fig. [6] are the spectra extracted from 90 per cent 
of the PSF, which corresponds to roughly a 50 arcsec width strip 
across the cluster. 



4.1 Spectral fitting 

The spectral fitting procedure here broadly follows that used by 
[Sanders et al.| ( |2008 ^. Given the wavelength range of the RGS spec- 
tra where the background in low, we allow the N, O, Ne, Mg, Ca, 
Fe and Ni metallicities to vary independently. We fix the metallic- 
ities of C, Si, S and Ar to the Fe metallicity, as we do not have 
good constraints for these elements. We use the PHABS absorption 
model, fixing the absorption column density to the Galactic value 
(as obtained from Kalberla et al.|2005[ shown in Table [TJ. We only 
fit the first order spectra in this analysis, using the wavelength range 
between 7 and 26A. The relative normalization of the RGS2 instru- 
ment was allowed to vary relative to RGSl. We grouped the spectra 
to have a minimum of 20 counts per spectral bin and minimized the 
in the spectral fit. We did not use the C statistic here because the 
template background spectra are synthetic count rate spectra, not 
count spectra, so Poisson statistics cannot be assumed. 

As clusters are extended objects, we need to account for the 
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Figure 4. Projected emission-weighted temperature maps of each object 
from Chandra data (right panels) next to X-ray images (left panel). The 
colour bars at the bottom of each pair of images show the temperature in 
keV. Note the ring of cool gas in Abell 262 to the north east (marked by 
arrows). Each region in the temperature map was chosen to have a signal to 
noise ratio of 22 (approximately 480 counts). 



spectral broadening due to the spatial extent. This is done using a 
smoothing model ( [Sanders et al.|2008"] ). We use an upper limit of 3 
arcmin for the broadening when fitting. Different broadening s are 
used for components where feasible. We caution that the metallic- 
ities from the spectral fitting may be low because as we assume a 
Gaussian broadening from the spatial extent and do not fully model 
the wings of the emission lines. Metallicity ratios are more likely 
to be accurate than the absolute values. 

It is possible to use the spatial information from imaging with 
Chandra to account for the spatial broadening of the emission lines, 
using for example the RGSXSRC model in XSPEC. This approach 
assumes that the spectral lines are emitted from the same region 
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Rest wavelength (A) 

Figure 5. (Top panel) The fluxed spectra as a function of rest wavelength. The stronger emission lines have been labelled. (Central panel) Model solar 
metallicity spectra at four different temperatures for comparison, broadened to comparable spectral resolution. (Bottom panel) Spectrum of a solar metallicity 
cooling flow model cooling from 2.5 to 0.0808 keV. 



as broad band imaging and that the X-ray spectrum is uniform over 
this region. This is likely to be an incorrect assumption if the cooler 
gas is more centrally concentrated as suggested by the temperature 
profiles and maps (Fig. [2] and Fig. |4]). We therefore use separate 
free parameters for the spatial broadening for the hotter and cooler 
components. Our spectral fitting of the RGS data confirms that the 
cooler gas is more spatially concentrated. The good agreement be- 
tween the simple spectral fitting used here and the smooth particle 
analysis in Section |5] which treats the spatial extent of the source 
self consistently, implies the simple broadening approach works 
well. 



where is the angular diameter distance to the object in cm, V is 
the volume in the cluster being examined in cm^ , z is the redshift 
and Hq and «h are the electron and Hydrogen number densities, 
respectively, in cm~^ . The units of are cm~^. 

The two temperature models are reasonable fits to the data, 
given the systematic issues with fitting high spectral resolution 
data. The lower temperature value for the three targets is close to 
0.7 keV. This component is also broadened less than the hotter com- 
ponents, being consistent with it being emitted from a smaller re- 
gion. 



4.1.1 Two temperature fits 

We first examined some fairly simple spectral models. As a sin- 
gle temperature model is clearly insufficient to fit the data we 
used a model consisting of two APEC components with the vari- 
able elemental abundances listed above (we refer to this model as 
2 X VAPEC). Each of the thermal components were allowed to have a 
different spatial broadening and the metallicities of the components 
were tied to be the same. The temperatures and model normaliza- 
tions were allowed to be free. The best fitting parameter values and 
uncertainties are shown in Table [4] The spectra with the best fitting 
models are shown in Fig. [T] with the residuals from the fit. 

In XSPEC the normalization of the thermal components (pro- 
portional to the emission measure) is defined as 

^ = ^ / ^e^HdV, (1) 



4.1.2 Power law emissivity model 

Another simple parametrization is to have a power law distribu- 
tion of emission measures as a function of temperature. The emis- 
sion measure of each temperature component is proportional to 
(r/Tmax)", where Tmax is a parameter for the upper temperature. 
This model is not physically motivated. We chose to have the spec- 
tra in the model constructed from APEC spectra and also included 
another APEC component to account for any non-cool gas. 

We show the best fitting parameters for this model (which we 
refer to as VAPEC-i-CEVMKL) in Table |4] In this model we find 
the best fitting indices for the power law distribution of 3.1 for 
Abell 3581 and HCG 62 and 2.6 for Abell 262. It is unclear whether 
this model is a useful parametrization as the relative balance be- 
tween the APEC and power-law emissivity models are quite differ- 
ent between the clusters. The WDEM power-law emissivity model 
( [Kaastra et al.||2004| with high and low temperature cut-offs ap- 
pears to be a good parametrization for many clusters, though we do 
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Table 4. Spectral fitting results. (1) Temperatures of the two components for 2xvAPEC, upper temperature of CEVMKL and VAPEC temperature for 
VAPEC+CEVMKL, Upper temperature of mkcflow and vapec temperature for vapec+mkcflow and lower temperature of mkcflow if there is a free 
minimum. (2) Emission measures of thermal or CEVMKL components multiplied by 10^. (3) Mass deposition rates in yr"^ for MKCFLOW components or 
(4) emissivity index a for CEVMKL components. (5) Smoothing scales in arcmin for first and second components, respectively. The N, O, Ne, Mg, Ca, Fe and 
Ni columns show metallicities in Solar units. Other metallicities are tied to Fe, except for He. 



Model 


Cluster 


Xv 


m (1) 
m (1) 


A (2) 
^2(2) 


M(3) 
a (4) 


sizei (5) 
size2 (5) 


N 
Mg 
Ni 


O 
Ca 


Ne 
Fe 


2 X VAPEC 


A262 
A3581 
HCG62 


1.20 = 
1488/1238 

1.14 = 
1674/1466 

1.18 = 
653.8/553 


0.72 ±0.02 
1.49 ±0.03 

0.69 ±0.03 
1.41 ±0.02 

0.69 ±0.01 
1.07 ±0.04 


6.9 ±0.4 
65.2±1.5 

4.3 ±0.5 

89±2 

10.0±0.9 
9.0 ±0.9 




0.34 ±0.06 
0.73 ±0.06 

Q 54+0.17 
0.71 ±0.04 

0.34 ±0.03 


0.23±0.17 
0.40 ±0.06 
1.47±0.15 
0.24±0.13 
0.43 ±0.05 
0.82±0.11 
0.54 ±0.21 
0.58 ±0.09 
0.8^^i 


0.25 ±0.02 
2.6 ±0.5 

0.24 ±0.01 
2.41 ±0.40 

n ooQ+0.036 
u.z,z,o_QQ2g 

0.7 ±0.5 


0.27 ±0.07 
0.37 ±0.03 

0.31 ±0.06 
0.32 ±0.02 

0.44 ±0.09 
0.43 ±0.04 


VAPEC+ 
CEVMKL 


A262 
A3581 
HCG62 


1.21 = 
1493/1238 

1.14 = 
1670/1466 

1.28 = 
705.8/553 


1.63 ±0.06 
1.49 ±0.04 
0.87 ±0.01 


37±4 
150 ±23 

65 ±1 

5.7 ±0.9 
81±9 


2.56±0.16 
3.1 ±0.5 
3.1±0.3 


1.7±0.3 
0.41 ±0.04 

0.9±0.2 
0.45 ±0.08 

>2.7 
0.33 ±0.03 


0.23±0.18 
0.44 ±0.08 
1.54±0.18 
0.26±0.14 
0.43 ±0.06 

r^ 77+0.13 

0.33±0.13 
0.56 ±0.08 
1.1±0.2 


0.28 ±0.03 
2.9 ±0.6 

0.25 ±0.02 
2.55 ±0.44 

0.15 ±0.02 
<0.56 


0.28 ±0.06 
0.42 ±0.04 

0.34 ±0.06 
0.35 ±0.02 

0.42 ±0.06 
0.33 ±0.03 


VAPEC+ 
MKCFLOW 


A262 
A3581 
HCG62 


1.24 = 
1534/1239 

1.15 = 
1683/1467 

1.35 = 
746.1/554 


1.44 ±0.03 
1.40 ±0.02 
0.83 ±0.02 


61±2 
87±2 
16±1 


4.4 ±0.3 
4.6 ±0.5 
3.6±0.3 


0.77 ±0.09 
0.43 ±0.06 

0.71 ±0.05 
56+0-^^ 

0.76±0.12 
0.19±0.05 


0.12±0.11 
0.41 ±0.06 
1.68±0.15 
0.18±0.11 
0.43 ±0.05 
88 + 12 
0.20 ±0.09 
0.50 ±0.07 
1.23±0.16 


0.19±0.01 
2.4 ±0.5 

0.21 ±0.01 
2.4 ±0.4 

0.10±0.01 
<0.27 


0.23 ±0.06 
0.32 ±0.03 

0.29 ±0.05 
0.31 ±0.02 

0.36 ±0.06 
0.29 ±0.02 


VAPEC+ 
MKCFLOW 
(free minimum) 


A262 
A3581 
HCG62 


1.20 = 
1489/1238 

1.14 = 
1672/1466 

1.22 = 
675.4/553 


1.59 ±0.05 
0.53 ±0.02 

1.45 ±0.03 
0.52 ±0.05 

0.96 ±0.02 
0.51 ±0.02 


53±2 
82±3 
6.7 ±1 


6.9 ±0.5 

8±1 
10.0±0.9 


^•^-0.1 

0.34 ±0.05 

0.74 ±0.06 

4Q+016 
^•^^-0.11 

>2.6 
0.34 ±0.03 


0.25 ±0.19 
0.41 ±0.07 
1.51±0.17 
0.25 ±0.14 
0.43 ±0.05 
0.80±0.12 
0.45 ±0.19 
0.57 ±0.07 


0.27 ±0.03 
2.8±0.6 

0.25 ±0.02 
2.5 ±0.4 

0.20 ±0.02 
<0.86 


0.28 ±0.07 
0.40 ±0.04 

0.32 ±0.06 
0.33 ±0.02 

0.42 ±0.06 

'^•-^'-0.02 


6 X VAPEC 


A262 
A3581 
HCG62 


1.19 = 
1464/1235 

1.10 = 
1602/1463 

1.39 = 
764/550 


See 

Fig.ji] 


See 

Fig.ji] 




0.81^^.08 
0.33 ±0.06 

0.72 ±0.05 
0.57±0.16 

1.7±0.7 
0.44 ±0.03 


0.29 ±0.19 
0.40 ±0.08 
1.59±0.19 
0.29±0.14 
0.40 ±0.06 
0.95±0.15 
0.42±0.18 
0.49 ±0.08 
1 04+0-1^ 

^•^^-0.23 


0.28 ±0.03 
2.9 ±0.6 

0.26 ±0.02 
2.4 ±0.5 

0.20 ±0.02 
<0.99 


0.32^0.10 
0.48^0-06 

0.34 ±0.07 
0.38 ±0.03 

0.41 ±0.07 
44+0-04 

^•^^-0.06 


6XVMCFL0W 
(APEC) 


A262 
A3581 
HCG62 


1.19 = 
1467/1235 

1.10 = 
1609/1463 

1.36 = 
747/550 


See 
Fig.[l4] 


See 
Fig.[l4] 




n 77+0.08 
'-0.06 

0.36 ±0.06 

0.73 ±0.05 
0.49±0.16 

> 1.9 
0.43 ±0.03 


0.31±0.19 
0.40 ±0.08 

1.69±0.19 

n OO+0.08 

'^•■^■^-0.14 

0.41 ±0.06 
0.88±0.15 
0.53 ±0.02 
0.62 ±0.09 

1 1+0.3 
^•^-0.2 


0.28 ±0.02 
3.0±0.6 

0.27 ±0.02 
2.8±0.5 

0.7 ±0.6 


31+0-09 
47+0-04 

'^•^'-0.02 

0.41 ±0.06 
0.44 ±0.02 

0.49 ±0.09 
0.48^0-03 
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Rest wavelength (A) 



Figure 6. Fluxed spectra, extracted from 90% extraction region and cor- 
rected to rest wavelength. 

not investigate that model here, because our later analysis shows 
that the temperature distribution in these cool objects does not ap- 
pear to be a powerlaw over a wide temperature range. 



4.1.3 Cooling flow models 

The spectrum can be fit by a cooling flow spectral model made up 
of a range thermal components where the emission measure distri- 
bution as a function of temperature is that of a gas which is radia- 
tively cooling in a steady state and with no heating. Here we use 
a VMKCFLOW model to calculate this. We use the APEC model to 
calculate its thermal components. 

We first fitted a model where the gas is cooling to the mini- 
mum temperature (0.0808 keV). The upper temperature, mass de- 
position rate and metallicities are free parameters. We also include 
the non cooling gas as an APEC component setting its temperature 
to the upper temperature of the VMKCFLOW model and using the 
same metallicities. We also allow for different spatial broadening 
for the thermal and cooling flow components. 

The best fitting parameters for this spectral model are shown in 
Table|4] hsted under VAPEC-hMKCFLOW. The mass deposition rates 
for the three objects are 4.4, 4.6 and 3.3 yr~^ for Abell 262, 
Abell 3581 and HCG 62, respectively, where the cooling flow min- 
imum temperature is fixed to zero. 

If we instead fit a cooling model where the minimum temper- 
ature is allowed to be non-zero, we obtain mass cooling rates of 
6.9, 9 and 10 yr~^ . The systems show very similar minimum 
temperatures of 0.53, 0.52 and 0.51 keV, for Abell 262, Abell 3581 
and HCG 62, respectively. 




Wavelength (A) 

Figure 7. Two component VAPEC fits to the RGS data extracted from each 
object. The models are the solid lines. The data are shown as points. The 
spectra and models from the two RGS detectors are plotted in black and ma- 
genta. Some spectral regions are covered by only one detector. The panels 
with the axis scale x show the residuals of the fits. The data and models are 
shown dividing through by the effective area of the instrument as a function 
of wavelength. The data have been rebinned to have a signal to noise ratio 
of at least 10 in each spectral bin. 

4.1.4 Multi-component thermal models 

To investigate more thoroughly the amount of gas as a function of 
temperature, we have fitted a model made up six thermal compo- 
nents at fixed temperature. We use temperatures of 0.272, 0.544, 
0.862, 1.366, 1.719 and 2.164 keV. These particular temperatures 
were chosen because they are temperatures the APEC model is cal- 
culated at in the table in XSPEC, and they are sufficiently sepa- 
rated in temperature in order to be resolved. The three coolest com- 
ponents were smoothed by a different spatial scale to the three 
hottest. We tried two different thermal spectral models, APEC and 
the MEKAL model taken from SPEX version 2.00.11. With this 
model we fix the temperature but allow the normalizations to vary 
to examine the emission measure distribution. 

Fig. |8] shows the normalization of each component as a func- 
tion of temperature for each of the clusters with the two the spectral 
models, APEC and SPEX. Also shown on the plot are the expected 
normalization distributions for a cooling flow of lOM© yr~^ , us- 
ing the same metallicities as obtained from the spectral fits. These 
lines vary because the objects lie at different distances and have 
different metallicities. At 0.862 keV, the objects have normaliza- 
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Figure 8. Normalizations for the temperature components in a 6 component spectral model. The left and right panels show the results for APEC and SPEX 
models, respectively. The solid lines show the expected distribution of normalization for a cooling flow of 10 M© yr~^ in each of the clusters. 



tions compatible with 10 yr~^ , but at lower temperatures the 
maximum radiative cooling rate is much lower. Note that the SPEX 
spectral model gives significantly lower normalizations of the 0.544 
and 0.272 keV components than APEC. 

The spatial broadening of the lines from the coolest gas imply 
they come from a regions of approximate radii 0.33, 0.57 and 0.44 
arcmin, for Abell 262, Abell 3581 and HCG 62, respectively. 

We can also compare these results against the normalization 
distribution we can measure from Chandra data. If we fit a multi- 
temperature model to individual spatial regions we can map the dis- 
tribution of gas as a function of temperature. We used the Contour 
Binning algorithm to choose spatial regions from the Chandra data. 
As the Abell 262 observation had more counts than the other two, 
we chose a minimum signal to noise ratio of 32 (~ 1024 counts). 
For the other two objects we used 22 (~ 480 counts). We extracted 
spectra from each of the regions. We fitted a model to each spec- 
trum made up of 0.544, 0.862, 1.366 and 2.164 keV APEC compo- 
nents. We did not include the 0.272 and 1.719 keV components as 
the data are not sensitive enough to measure these. The metallicity 
was allowed to vary in each fit, but the metallicities of the tempera- 
ture components were tied together. The absorbing column density 
was tied to the Galactic value. The C statistic was minimized in the 
spectral fitting. 

Fig.[9]shows the resulting normalization per unit areas for each 
of the objects in each of the temperature components. The images 
show we detect cool temperature components in the core of each 
cluster as with the RGS data. 

The volume filling fraction of the coolest temperature compo- 
nent can be calculated assuming that the components are in pres- 
sure equilibrium. The volume filling fraction of component / is 



where is the normalization of component /, 7] is its temperature 
and j sums over all the components. We show the volume filling 
fraction for the coolest component in Fig.^j This fraction is very 
low for Abell 262 and Abell 3581, indicating this is true multiphase 
gas. The value in HCG 62 is much larger. The projected tempera- 
ture in the core of HCG 62 lies between 0.544 and 0.862 keV so the 



signal in these components there is probably not due to multiphase 
gas. There could still be some multiphase gas at lower temperatures 
still. 

We also calculate the mass in the 0.544 keV components as- 
suming they are in pressure equilibrium with the pressure of the 
gas calculated using the Chandra deprojection analysis. These val- 
ues are plotted on Fig|3] compared with the distributions from the 
Chandra radial analysis. 



4.1.5 Markov Chain Monte Carlo spectral analysis 

The fixed temperature component models provides a reasonable de- 
scription of the data. To construct a normalization plot without dis- 
crete temperature bins we used a MCMC analysis. 

We constructed a model made up of four APEC components 
with the metallicities fixed to be the same. These components were 
first fixed to be at 0.5, 0.8, 1.2 and 2 keV. The coolest three com- 
ponents shared the same spatial broadening size. The 2 keV com- 
ponent was broadened by a separate value. We fit the model to get 
initial metallicities, normalizations and spatial broadening sizes. 

In the MCMC analysis, we included in the chain the tem- 
peratures of the three lowest temperature components, their single 
broadening size and the metallicities of all the components. To con- 
struct the chain we used a custom proposal distribution in XSPEC 
which dynamically adjusted a Gaussian proposal distribution for 
each parameter in every 100 iterations so as to have a repeat frac- 
tion close to 75 per cent during the burn in period. Each parame- 
ter was stepped sequentially in the chain in random order. We ran 
the chain three times for each object. Each chain had a length of 
10^ and the initial 2x10^ values were discarded. The temperatures 
were constrained to lie between 0.08 and 4 keV. 

Fig. [TT] shows the distribution of temperatures and emis- 
sion measures for the three unfixed temperature components. The 
fraction shown on the image is the fraction of times one of the 
free temperature components in the chain was in that particular 
temperature-normalization bin. There are 40 logarithmic temper- 
ature bins between 0.09 and 2.9 keV and 40 logarithmic normal- 
ization bins between 10~^ and 10~^ cm~^ plotted. We also include 
the results from running the MCMC analysis on a simulated 120 ks 
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0.544 keV Abell 262 

20 kpc (60 arcsec) 
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0.544 keV Abell 3581 

20 kpc (45 arcsec) 




1.366 keV 




Figure 9. Normalization distributions as a function of temperature from the Chandra imaging. The colour scales are logio XSPEC normalizations per square 
arcsec. The regions for Abell 262, Abell 3581 and HCG 62 contain a signal to noise ratio of 32, 22 and 22, respectively. 
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Figure 10. Maps of the volume filling fraction of the 0.544 keV component in the fits to Chandra data, calculated from the values in Fig. [9] 



dataset for a lOM© yr~^ cooling flow cooling from 2 keV at the 
redshift of Abell 262 and with 0.6 Z© metallicity. The dataset was 
simulated with the Abell 262 RGS responses and a single spatial 
broadening. 

These distributions are not completely straightforward to in- 
terpret. The results from the simulated dataset show that fitting dis- 
crete temperature components to a spectrum generated from contin- 
uous distribution yields particular temperatures, rather than a con- 



tinuous distribution. The gas temperatures which are easily differ- 
entiated spectrally in the RGS data are likely to be these particular 
temperatures. 

In the real data, the high likelihood temperature-normalization 
regions have similar temperatures to the two hottest likely regions 
from the simulated data. The coolest region in the simulated data 
(around 0.3 keV) is not seen in the real data. In the real data the 
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Figure 11. MCMC distributions of temperature and normalization from 
a multitemperature model with three variable temperature components 
and one fixed at 2 keV. Also shown is the distribution from a simulated 
lOM© yr~^ cooling flow at the redshift of Abell 262. The points plotted 
are the results from the 6 component VAPEC fits in Fig. [8] 



third temperature component moves around in the allowed temper- 
ature space giving the continuous distribution seen. 

The simulated example shows that any "special" temperatures 
that come out of an analysis should be treated with caution. A spec- 
trum generated from a continuous distribution of temperatures will 
yield particular temperatures under analysis as a matter of course. 
These particular temperatures will depend on a combination of in- 
strumental response and how spectral features change as a function 
of temperatures. 

The lack of the third component, seen in the cooling flow 
model simulation results around 0.2 keV, in the real data confirms 
that we do not observe any gas around these temperatures. We also 
see very good agreement between the simple 6 component spectral 
fitting results from S ection [4 . 1 . 4 | (Fig . [8|l which are also plotted as 
points in Fig.[TT] and these MCMC results. 

If we increase the number of temperature components from 4 
to 6 in the MCMC analysis, we produce plots of temperature and 
emission measure shown in Fig. [12] The results are very similar to 
the previous plot for real data, showing the results are robust. For 
the simulated data we see a continuous distribution between the 
three emission measure peaks. 



4.1.6 Redshift from spectral fitting 

We can also measure the mean redshift of the intracluster medium 
in each object by spectral fitting. Using the 6x VAPEC model from 
Section [4. 1.4| we allowed the joint redshift of the 6 thermal com- 
ponents to be free. 

» shows the change in the quality of the fit {x^) as a func- 



Fig, 



13 



tion of redshift. We find the redshifts we obtain are roughly consis- 
tent with the measured optical values in Abell 3581 and HCG 62, 
but not in Abell 262. There is an offset of around 360kms~^ be- 
tween the two values. Such a shift could be due to in incorrect posi- 



Figure 12. A 6 component MCMC distribution of temperature and normal- 
ization. This can be compared against the 4 component results shown in 

Fig.[TT] 



tion for Abell 262 during the processing. However the wavelength 
shift at 19A would require a large offset of around 8 arcsec. There- 
fore it is probably a real offset between the optical measurement 
and X-ray value. [ S^ai et al.| ( [l994| ) examine the velocity distribu- 
tion of the galaxies in Abell 262 as a function of the distance from 
the cluster centre. Our velocity is within their distribution of ve- 
locities, but not at the peak. There is no trend towards our value 
towards the cluster centre. The D galaxy NGC 708 lies at the peak 
of the X-ray emission ( |Jones & Forman|1984| ), but its velocity does 
not match our velocity measurement either. There may be a bulk 
flow of gas in the cluster giving our velocity offset. 

We did not use the best fitting redshift of Abell 262 when ex- 
amining the temperature distribution of the ICM. If we allow the 
redshift to be free when fitting the six thermal components spec- 
tral model (Section fi.l.4| ), it makes no significant changes to the 
normalizations of each component. The error bars become slightly 
larger but there is no systematic shift in the best fitting values. 



4.1.7 Multi-component cooling flow spectral fitting 



As an alternative to the multitemperature model from Section [4!l.4| 
we can fit models which parametrize the amount of gas which could 
be cooling as a function of temperature. We used a model made up 
of several components which modelled gas cooling through differ- 
ent temperature ranges. We firstly used six bins between tempera- 
tures of 3.88 to 1.94, 1.94 to 1.37, 1.37 to 0.97, 0.97 to 0.49, 0.49 
to 0.24 and 0.24 to 0.080 keV. 

Fig. [14] compares the mass cooling rates obtained using the 
APEC spectral model and the SPEX spectral model. The results are 
fairly similar between the two models except in the 0.49 to 0.24 
bin and for the different limits in the coolest temperature bin. Also 
shown on the plot are the cumulative mass deposition rates as a 
function of temperature obtained with the Chandra surface bright- 
ness deprojection (Fig.|2]). The cumulative mass deposition rate at 
the cooling radius is shown as a shaded bar (the values are taken 
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Figure 13. Change in quality of fit as a function of redshift for each of the clusters from the RGS spectra. These results were from the 6 component APEC 
model with free redshift. The vertical line shows the redshifts taken from Table^ The horizontal lines show the A^^ values appropriate for confidence regions 
of 1, 2 and 3cr. 
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Figure 14. Comparison of mass cooling rates obtained with APEC and 
SPEX cooling flow models using 6 temperature bins. The shaded regions 
show mass deposition rates derived from Chandra surface brightness de- 
projection (Table [Sj. The solid line was created by plotting temperatures 
against cumulative mass deposition rates from Fig.[2| with the temperatures 
smoothed with a 3 bin sliding average. 



Figure 15. Comparison of mass cooling rates obtained with APEC cooling 
flow models using 6 or 9 temperature bins. The shaded regions show mass 
deposition rates derived from Chandra surface brightness deprojection (Ta- 
ble [Sj. The solid line shows temperatures plotted against cumulative mass 
deposition rates from Fig. [2] 



from Table [3]). There is quite a good agreement between the shape 
of the RGS temperature-mass cooling rate relation and the Chandra 
profile version, except for an offset in temperature between them. 



Fig. [15] shows a comparison between the results using the 
APEC model with 6 bins and one with 9 temperature bins. The 
results are consistent between the two different temperature bin- 
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Figure 16. RGS spectra around the expected position of the Ovii lines. 
These spectra have not been fluxed but have been background subtracted. 
The thick solid line is a model cooling flow spectrum cooling from 3 keV 
to zero temperature at the rate of 10 M© yr~^ (at Solar metallicity). 

nings. As expected, the constraints are relaxed with more tempera- 
ture bins. 

In Fig. [16] we show unfluxed background subtracted spec- 
tra between 21 and 23. 5 A, showing the region where O Vll lines 
should be observed if they exist. Also shown are model spectra for 
a 10 M0 yr~^ cooling flow for Solar metallicity material. O VII is 
one of the strongest indicators of gas around 0.25 keV in our wave- 
length range, however it is not one of the main coolants. Therefore 
the strength of the line will depend on the metallicity of oxygen and 
amount of cool gas. 

We can examine how the constraint on the upper limit on the 
amount of gas cooling in the lower temperature bin is affected by 
the oxygen metallicity. Lower oxygen abundances could potentially 
be caused by a relative deficit of Type II supemovae in the core of 
the clusters. In Fig.[T7]we show confidence contours for the mass 
cooling rate as a function of oxygen metallicity, for the lowest tem- 
perature bin in the APEC cooling model (see the best fitting values 
in Fig. [14]). The plots show that we could have larger amounts of 
cooling at low temperatures if oxygen metallicities are depleted at 
these temperatures, but very low metallicities are required to match 
the cooling rates seen at higher temperatures. 

Fig. [Ts] shows a comparison of the HCG 62 data with a model 
cooling to the minimum temperature at 6.7M0yr~^ at 0.5 Z©. 
This model is broadly consistent with the Fe XVII emission lines 
close to 15 and 17 A, but is inconsistent with the O VIII emission 
line at 19A and the O VII lines at 22.6 and 22.1 A. Reducing the 
oxygen abundance would help the fit, but we would need unrealis- 
tically low values (~ 0.05 Zq) to match the data. 

4.2 Line ratios 

A different way of examining the temperature distribution is to look 
at temperature sensitive emission lines, such as the different Fe 
XVII emission lines. In our spectra the strongest lines are the 3C 
line at 15.01 A (a 2pi^2 ~ '^^3/2 transition) and the combined 3G 
and M2 lines at 17.1 A (both 2p — 3s transitions). The 3F line at 
16. 8 A (another 2p — 3s transition) is also clearly seen in HCG 62. 
The 3D line at 15.26A is blurred into the 3C lines. The useful 
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Figure 17. Confidence contours on limits mass deposition rate in the 0.24 
to 0.08 keV temperature range as a function of oxygen metallicity. This was 
measured with the APEC model (see Fig.[T4}. 



3C/3D ratio cannot be used because of the blurring and low signal 
to noise of the 3D line. The 3C line can also be subject to resonant 
scattering due to its oscillator strength. 

We can measure the 3C/(3G-i-M2) ratio. Unfortunately there 
are disagreements between theory and experiment over how this 
ratio varies as a function of temperature (Bei ersdorfer et al.|2002| 
|2004| ). In addition, the cluster is likely to contain a range of gas 
temperature, so the line ratio will not come from a single temper- 
ature component. Nevertheless, we measured the strength of the 
emission lines. 

These results were measured by first taking the 6 component 
APEC spectral fit from Section |4.1| We then modified the APEC 
table of lines at each temperature to make the emissivity of the 
Fe XVII and O Vll lines to be zero. We then added four redshifted 
Gaussian components for the Fe XVII lines and two for the O Vll 
lines. The width of the Fe XVII lines were constrained to be the 
same. The width of the O Vll lines were set to be zero. We fit the 
spectrum over the same range of wavelength as in Section [4T7] The 
line strengths or upper limits we obtained are shown in Table [5] 
This spectral fitting procedure was used to make sure that we had a 
good model for the continuum around the lines. A simple powerlaw 
model was not sufficient, giving results dependent on the range of 
wavelength fitted. 

The 3C/(3G-hM2) ratio for each cluster is plotted in Fig. [19] 
In addition we plot the line ratio as a function of temperature from 
APEC, SPEX, MAPEC ( |Gu]|2007| ) and CHIANTI ( |Dere et al.]|1997[ 
|Landi et "aL]|2006| ). Most of the models show similar ratios as a 
function of temperature, except for APEC. 

|Doron &^ehar ( 2002 ) presented a detailed study of the Fe 
XVII emission line ratios. They included the effects of radiative re- 
combination, dielectronic recombination, resonant excitation and 
inner- shell collisional ionization, in addition to the collisional exci- 
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Figure 18. HCG 62 spectrum and comparison with a 6.7 yr cooling flow (cooling from 1.2 to 0.0808 keV with 0.5 Z© metallicity). 



Table 5. Measured line fluxes (10 erg cm ^ s ^ ). These are measured from a 95 per cent extraction region. Uncertainties and limits are la. 
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Figure 19. Measured Fexvii line ratios (in terms of energy flux) compared 
to theoretical models. The ratio measured in a cluster is shown by a line 
surrounded by a shaded region which represents the 1 a uncertainty. 



tation. They found that the 2p-3s lines at 16.78, 17.05 and 17.10A 
were enhanced by 25, 30 and 55 per cent, respectively. APEC is the 
outlying spectral model in Fig. [19] and does not include the results 
of |Doron & Behar|p002| ), though work is being done to improve 



the Fe XVII line strengths (R. Smith, private communication). SPEX 
does include the results of |Doron & Behar| ( |2002| ). The line ratios 
are not strong functions of temperature above ~ 0.5 keV, so it is 
difficult to get an exact temperature measurement using Fexvil. 

The APEC model, generally produces a better quality fit over- 
all compared to SPEX. The two component thermal fit (2xVAPEC) 
for Abell 262 has a lower by 100 compared to the SPEX model. 
We therefore use APEC for most of the analysis in this paper. 

We do not see any evidence for resonant scattering. However, 
we would not expect to see this necessarily, as the region where 
any scattering would take place is likely to be inside our spectral 
extraction region. The total spectrum from this region would re- 
main the same if there is only scattering within it. Analysis of the 
line profiles is required to look for scattering (e.g. |Xu et al.|2002| . 



5 SMOOTH PARTICLE INFERENCE MODELLING 

Rather than the traditional analysis we have performed in Section|4] 
we have also applied a smooth particle inference analysis method 
to the XMM-Newton data, as described in Peterson et al.| ( |2007| ). 

This smooth particle analysis is more complete than the 
MCMC procedure we outlined in Section |4.1.5| In that analysis 
we attempted to model the RGS spectrum with just the sum of 
4 or 6 spectral components. Here the whole intracluster medium 
is modelled as a set of X-ray emitting smoothed particles. Each 
particle is described by a set of parameters including temperature, 
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metal abundance, size and location. Also included in the model are 
components which account for the instrumental and cosmic X-ray 
backgrounds. The X-ray photons from the model are then propa- 
gated through the response of the X-ray detector and compared to 
the real data using a MCMC. Once the chain has converged prop- 
erties about the cluster can be extracted from the chain. 

This method is powerful as it can be applied to data from 
different kinds of instruments, such as gratings and CCD imaging 
spectroscopy. It builds up a model which fits the spatial and spectral 
data simultaneously. 

We use particles which are modelled as Gaussians spatially, 
with parameters for the temperature, metal abundance, redshift, 
normalization, absorbing column density and position and width 
on the sky. The abundances are restricted to lie between and 
2Z0. The Gaussian widths of the particles can lie between 0.03 
and 4.08 arcmin. The prior for the distribution of particle temper- 
atures is constant in — 1 < Xogi^ikT /keM) < 1. The redshift of the 
particles and their absorbing column density are fixed. An isother- 
mal absorbed MEKAL model is used to generate the X-ray spectrum 
from each particle. 

The soft X-ray background is modelled as a spatially uniform 
unabsorbed MEKAL component fixed at 0.15 keV (for Abell 262 
and Abell 3581) and 0.3 keV (for HCG 62). The component is fixed 
at z = and Solar metallicity. The harder X-ray background, pro- 
duced primarily by unresolved extragalactic sources, is modelled 
as a r = 1 .47 spatially-uniform powerlaw, with Galactic absorp- 
tion. The detector background is also modelled by emission lines 
representing the contribution from high energy particles and elec- 
tronic noise, modelled as an exponential. The detector noise is fixed 
between the different observations. 

The photons emitted by the particle and background models 
are propagated through the RGS and EPIC responses and compared 
with the input data. The fraction of photons emitted by each compo- 
nent and the relative normalizations of the models are free param- 
eters. Energy cuts of 0.3 to 10 (EPIC MOS), 1.1 to 10 (EPIC PN) 
and 0.2 to 2.5 keV (RGS) are apphed to the photons. The MCMC 
is run for 1 000 (EPIC) or 2 500 (RGS) iterations. The output of 
every iteration in the chain after 200 is used to produce the set of 
model parameters. For this analysis we examine only the longer 
XMM-Newton datasets for each object (observations 0504780101, 
050478301 and 0504780501), except for the EPIC observation of 
Abell 262, where we use 0504780201 due to the flaring in the 
longer observation 0504780101. 

In Fig.[20|are shown the resultant distribution of particle tem- 
peratures from the MCMC chains. The total normalization of the 
particles in each logarithmic temperature bin is plotted divided by 
the bin width. Also plotted for comparison of the 6 component 
SPEX results from Fig. [8] divided by an approximate bin width 
(there are no exact bin boundaries between the temperature compo- 
nents or at the edges of the distribution). We plot the SPEX spectral 
fitting results as the spectral model is very similar to the MEKAL 
model used in the smoothed particle inference analysis. 

Below 2 keV temperature, there is very agreement between the 
RGS smooth particle results and the simple SPEX spectral fitting. 
At high temperatures the results disagree as the simple spectral fit- 
ting only includes components up to 2.164 keV in temperature and 
examines the RGS data from only the central region of the cluster. 
The EPIC instruments are more sensitive to hot gas than the RGS 
detectors and so it is not surprising these results disagree with the 
RGS smooth particle values at high temperatures. Below 0.5 keV 
there is little sensitivity to cold gas in the EPIC detectors. 




0.1 



Temperature (keV) 

Figure 20. Results from a smooth particle modelling of the clusters, shown 
as solid (RGS) and dotted (EPIC) lines. Also plotted for comparison are 
the simple 6 component SPEX spectral fitting results from Fig. [S] (plotted as 
points). Emission measures have been divided by temperature bin widths. 



6 DISCUSSION 

6.1 The nature of the coolest X-ray emitting gas 

Purely from the existence of Fe XVII emission lines in our spectra 
of these objects we can see that they contain X-ray emitting gas 
somewhere between approximately 0.15 and 0.80 keV (Fig. |2TJ. 
Using spectral fitting, the APEC and SPEX models indicate that there 
is significant gas at around 0.544 keV (Fig. [8]). 

We can calculate the implied cooling time of the coolest com- 
ponent using its measured emission measure. We assume that the 
cool gas is in pressure equilibrium with its surroundings, using the 
peak thermal pressure from the Chandra profiles (Fig. [5]). Taking 
this pressure and temperature we compute the density and with the 
fitted metallicities we estimate the cooling time. The cooling times 
are shown in Fig. [2] From the density and the emission measure we 
can also estimate the volume of the coolest 0.544 keV component. 
These volumes are around 74, 380 and 612 kpc^ for Abell 262, 
Abell 3581 and HCG 62, respectively. Using the width of the spec- 
tral lines we can also estimate the area on the sky from which they 
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Figure 21. Ionization fraction of different iron ions as a function of temper- 
ature. Data are taken from |Mazzotta et aL]jl998) . 

are emitted (Section [4T| . If the gas is volume filling, the volume 
should be compatible with the area. However, from the line widths, 
the implied volume filling fractions for the 0.544 keV component 
of 6, 3 and 31 per cent, for Abell 262, Abell 3581 and HCG 62, 
respectively. 

Single phase analyses of CCD data miss this very coldest gas. 
An intermediate temperature between the volume filling compo- 
nent and the cold blobs is obtained. This can be seen in the discrep- 
ancy between the Chandra profiles and the RGS values in Figures 
[T4]and[T5] In the Perseus cluster multiphase analysis is required to 
detect the filaments spectrally (figure 12 in Fabia n et al.|2006| ) and 
in 2A 0335-1-096, where the multiphase analysis was required to 
reveal the morphology of the coldest X-ray emitting gas ( [Sanders | 
let al.|2009a| ). 

The Chandra volume filling fraction maps (Fig. [T0| of the 
coolest components give very similar results to the calculations 
from the RGS results. The small volume filling factor and spatial 
extent (from Chandra maps and from the line widths in the RGS 
data) for Abell 262 and Abell 3581 (Fig. [9} suggests that these 
components consist of small cool blobs in a hotter medium. It is 
a true multiphase gas. In HCG 62 the volume filling factor of the 
0.544 keV component is much larger (around 30 per cent). This 
temperature is also close to the mean projected temperature from 
the cluster, so we are probably detecting the volume filling phase 
(which lies in temperature between 0.544 and 0.862 keV). The cool 
blobs in Abell 262 and Abell 3581 are likely to be similar to those 
found in other clusters, such as 2A 0335-1-096. 

The temperature distribution in HCG 62 is also narrower than 
in the other two objects, as seen from spectral fitting (Figures [8] 
and[T4|, the simple MCMC analysis (Fig.[TT} and the full smooth 
particle inference analysis (Fig.[20|. 

6.2 The distribution of gas temperature 

[Peterson et al.| ( [20 0T) found that the differential luminosity distri- 
bution of gas as a function of temperature in a sample of clusters 
showed a powerlaw index of 1-2 up to an upper temperature. To 
examine whether this holds in these objects we show the distribu- 
tion of gas as a function of temperature relative to a cooling flow 
model in Fig. [22] To calculate this, we took the gas normalization- 
temperature distribution from the smooth particle analysis (Fig.[20| 
and divided the normalizations by those from a cooling flow model 
at each cluster redshift. The cooling flow model was calculated 
with Solar metallicity and with the mass deposition rate taken from 
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Figure 22. Distribution of gas relative to the amount expected from a cool- 
ing flow model. The cooling flow mass deposition rates used were taken 
from the surface brightness deprojection results at the cooling radius (listed 
in Table [3}. The top panel shows the gas distribution as a function of tem- 
perature relative to the temperature at the cooling radius. Powerlaw distribu- 
tions of indices 1 and 2 are shown for comparison. The bottom panel shows 
the distribution as a function of absolute temperature. 

the Chandra surface brightness deprojection at the cooling radius 
(listed in Table [3]). The top panel in Fig. [22] shows the distribution 
as a function of temperature relative to the temperature at the cool- 
ing radius. The bottom panel shows the distribution as a function 
of absolute temperature. 

In these objects the distribution appears more complex than a 
simple powerlaw distribution with an index of between 1 and 2. If 
plotted at absolute temperature there appears to be a sharp break at 
around 0.5 keV, with much less gas below this temperature. This 
break is consistent with the results from a cooling flow model with 
a non-zero minimum temperature (Section [4. 1.3| l. A powerlaw is 
consistent with the results above this temperature. The results of 
the cluster sample in Peterson et al.| ( [2003 ) have no detections of 
gas below 0.5 keV (in figure 6 in that paper). These results would 
also be consistent with a powerlaw distribution above 0.5 keV and 
a break below that temperature, although the data quality for many 
of the objects makes this hypothesis difficult to test here. A tem- 
perature break around 0.5 keV was also seen in M87 ( [Werner et aL[ 
[20061 [Simionescu et al.|2()08] ). 
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If the breaks in the distributions around 0.5 keV are 
widespread, this could indicate that an extra physical process be- 
comes relevant below around 0.5 keV temperature. This could, for 
example, be the onset of non-radiative cooling (see Section [64] be- 
low). Such a break would be easier to see in cooler objects such as 
we are examining in this paper. At this threshold temperature the 
X-ray emitting gas has cooling times less than 10^ yr. 



6.3 Heating the coolest gas 

All three of these clusters show signs of AGN activity in their cores. 
If AGN activity is supplying the correct amount of heat to prevent 
cooling, the problem is how can the AGN provide the right amount 
of heat to the coolest X-ray gas, which only occupies a small frac- 
tion of the core of the cluster? This material has the shortest X-ray 
cooling time and lowest temperature. Any sort of heating mecha- 
nism which heats in a volume- averaged way should not work. The 
heat would need to be targeted to this X-ray coolest gas in order to 
prevent that cooling. 

The cool blobs are around half the temperature of their sur- 
roundings so their emissivity is about four times larger (the cool- 
ing curve is approximately flat over this temperature range). They 
therefore require a heating rate four times larger than their sur- 
rounds. If there is widespread heating per unit volume then the sur- 
roundings would be overheated to prevent the coolest blobs from 
cooling further still. 

The material we observe in these cool clusters and groups is 
only a factor of 2.6 to 5 times smaller than the gas further out. The 
heating tuning problem is even worse for other objects. As men- 
tioned in the introduction, there is cool gas a factor of eight smaller 
than the outskirts in 2A 0335-^096 ( [Sanders et aTr 2009a), fifteen in 
Abell 2204 { [Sanders et al.|2009b| ), ten in Centaurus (Sand ers et al.j 
[2008] ), and twelve in the filaments in Perseus ( [Fabian et al.|2006j ). 



6.4 An alternative to heating: non-radiative cooling 

If AGNs are responsible for the bulk of the heating in clusters, then 
either the cool gas must be heated by other means, a mechanism 
is in place for channeling AGN heating to this cool gas, or there 
are significant amounts cooling in the central regions. Significant 
amounts of cooling are possible if the gas is not cooling radiatively. 
For instance, mixing of X-ray emitting material with cooler gas 
would allow the material to cool much quicker ( [Fabian et al.|2002| ). 
Energetically, the amount of energy lost in the cooling could easily 
be lost in the infrared output of a central cluster galaxy. 

The missing luminosity in X-rays for lOM0yr~^ cooling 
from 1 keV at 0.6 Z© is around 2 x 10"^^ erg s~^ Abell 262 has an 
infrared luminosity of 8 x 10^^ erg s~^ (O'Dea et al. 2008 ), com- 



patible with the energy which would need to be emitted if cooling 
were taking place by mixing. 

Infrared emission is also seen from HCG 62 ( [Johnson et al.| 



|2007 ') with a peak lAjim flux of around 10 mJy from the galaxy 
at the X-ray centroid, NGC 4778. If the 8 and 24 jim fluxes are 
linearly interpolated to get the flux at 15 jim, the infrared luminos- 
ity of the galaxy can be estimated to be 1.4 x 10^^ erg s~^ , using 
equation 13 of [Elbaz et al.[ {2002). This is again, much more than 
the energy that would need to be emitted if the X-ray gas was cool- 
ing by mixing. 

In objects such as Perseus and Centaurus, cool X-ray emitting 
filaments are associated with material at still colder temperatures 
and Hoc nebulosity (e.g. [Fabian et al.|2003[ [Crawford et al.|2005[ 



[Salome et aL|[2006[ ). In 2A 0335-h096 Hoc emission is also seen 
associated with the X-ray cool blobs ( [Sanders et al.[[2009a| . The 
cool X-ray emission and colder material must be intimately related. 



6.5 Cooler gas and star formation 

The infrared luminosity in Abell 262 can be con verted to a star for- 
-1 ' ' ' 



mation rate of 0.5 yr~^ ( [O'Dea et 1^20081 ). If these stars were 
forming constantly at steady state from cooling X-ray gas, it is con- 
sistent with our lowest mass deposition rate upper limits. The clus- 
ter also has an Hoc emission with luminosity of 1.1 x 10^^ erg s~^ 
|Hatch et al.|2007t . In addition 4 x 10^ of H2 gas was detected 
in the central cluster galaxy ( [Prandoni et al.[2007 ). Both the ionised 
and molecular gas have similar velocity structures, appearing to be 
emitted from a rotating disc. 

Abell 3581 has extended Hoc emission ( [Danziger & Focardi 
1988 ). HCG 62 shows very weak Hoc emission ( [Spavone et al. 
2006 ). [Variuri & Anupama[ ( [199 6') calculate that the Ha emission 
relative to the continuum is smaller than for other traditional cool- 
ing flow objects. 



6.6 Mass deposition rates 

Like the previous analyses of RGS spectra of cool core clusters, 
we find significantly less gas at low temperatures than would be 
expected if cooling were taken place at a steady state in the absence 
of heating. The best example of this is HCG 62, where our limits in 
the 0.24 to 0.49 keV temperature range are two orders of magnitude 
lower than the cooling rates compatible with the surface brightness 
profile or spectra of the hot gas. In Abell 262 and Abell 3581 the 
spectra are consistent with higher fractional rates of cooling. 

The picture is made more complex by the discrepancy be- 
tween the different plasma models over the strength of the Fe XVII 
emission lines (Fig. [T9| and generally in spectral fitting (Figures 
|8][T4]and[T5]). In addition the iron and oxygen metallicities have 
some effect (Fig.[T7j figure 17 in [Sanders et al. 2008). If, for some 
reason, the coolest gas is metal poor, the cooling rate may be being 
underestimated. Unless they are very underabundant in metals, this 
probably could not account for the very low cooling rate seen in 
HCG 62. 



7 CONCLUSIONS 

The deep XMM-Newton RGS spectra of these objects show Fe XVII 
line emission, indicating material around 0.5 keV in temperature. 
The width of the emission lines and the spatial extent of the cold 
gas seen by Chandra show that the regions they are emitted from 
are resolved. Assuming pressure equilibrium, we can see that this 
cool component is not volume filling in Abell 262 and Abell 3581 
and likely consists of cool blobs such as those seen in 2A 0335-1-096 
or filaments as in Perseus or Centaurus. 

The strength of these Fe XVII emission lines and lack of O Vll 
emission lines show that there is substantially less material seen 
at lower temperatures than expected in the cooling flow picture, as 
seen in previous studies. 

The multiphase nature of this cold material presents prob- 
lems for heating mechanisms which are preventing the bulk of the 
material from cooling. The coolest material requires four times 
more heat than its immediate surroundings. Any heating mecha- 
nism stopping this material from cooling needs to target it in some 
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way, or it is cooling, perhaps non-radiatively. Non-radiative cool- 
ing, perhaps by mixing, could be important at the lowest tempera- 
tures. In Abell 262 and HCG 62 the infrared fluxes of the central 
galaxy could contain the heat output of this process. 
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